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ABSTRACT 

Context. The compact radio source Sagittarius A* (Sgr A*) at the centre of our Galaxy harbours a supermassive black 
hole, whose mass (~ 3.7 x 10 6 Mq) has been measured from stellar orbital motions. Sgr A* is therefore the nearest 
laboratory where super-massive black hole astrophysics can be tested, and the environment of black holes can be 
investigated. Since it is not an active galactic nucleus, it also offers the possibility of observing the capture of small 
objects that may orbit the central black hole. 

Aims. We study the effects of the strong gravitational field of the black hole on small objects, such as a comet or an 
asteroid. We also explore the idea that the flares detected in Sgr A* might be produced by the final accretion of single, 
dense objects with mass of the order of 10 20 g, and that their timing is not a characteristic of the sources, but rather 
of the space-time of the central galactic black hole in which they are moving. 

Methods. The problem of tidal disruption of small objects by a black hole is studied numerically, using ray-tracing 
techniques, in a Schwarzschild background. 

Results. We find that tidal effects are strong enough to melt sufficiently massive, solid objects, and present calculations 
of the temporal evolution of the light curve of infalling objects as a function of various parameters. Our modelling of 
tidal disruption suggests that during tidal squeezing, the conditions for synchrotron radiation can be met. We show 
that the light curve of a flare can be deduced from dynamical properties of geodesic orbits around black holes and that 
it depends only weakly on the physical properties of the source. 
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1. Introduction 

The centre of our Galaxy harbours the nearest 
massive bl a ck h ole. Stellar orbits determinations 
(|Ghez et al.l I2005T ) revealed a central dark mass of 
(3.7 ± 0.2) x 10 6 [i? /(8kpc)] 3 M©, with an updated value 
of (3.61 ± 0.32) x 10 6 M (7) from observations with SINFONI 
l|Eisenhauer et al.ll2005D confined within a radius of 45 AU. 
The proximity of the Galactic centre (~ 8 kpc) allows us 
to stud y the environmen t of massive black holes in detail 
(see e.g. lAlexanderll2005l) . 

Although the black hole at the Galactic centre is not 
an active one, flares from its direction have been detected 
in X-rays (Chandra and XMM-Newton satellites) and 
infra-red (VLT adaptive optics imager NACO, SINFONI 
infra-red adaptive-optics integral spectrometer on the ESO 
VLT). These flares are modulated on a sh ort timescale 
with average periods of ~ 20 minutes (iGenzel et al 



i ng ho t blobs (e.g. lYusef-Zadeh et alj (l2006D.lMarrone et alj 
(200 81)), hot spot/ ri ng models (e . g. iTrippe et alj (120071), 
lEckart et all (I2004D. feckart et all (l2006al). IBelanger et all 
(I2006D. iMever et all (l2006ah. iMever et all (l2006bD . 



Broderick fe Loebl (|2006h . IZamaninasab et alj (1 200 



20031: IBelanger et al.ll2006l : lEckart et alj2006bHMever et al 
2006a). Multi- wavelength campaigns found that the time 



lag between X-ray and NIR flare emissi on is small, strongly 
suggesting a common physical origin (lEckart et al.l l2004 
l2006atlYusef-Zadeh et al.ll2006t pTYippe et al.ll2007f) . 

Several different models were p roposed to expla i n thes e 
phenomena: d isk instabilitie s (e.g. jTagger fc MeHal (2006), 



lYuanl (120061) . iFalanga et all (120081). lEckart et all (120081)) 



star-disk interaction (e.g. iNavakshin et alj ( 20041 )). expand 
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The article by I Alexander] (120051) contains an exhaustive 
description of the Galactic centre (GC) environment and 
its stellar dynamics. It is reasonable to expect that stars at 
the GC are surrounded by planets and low-mass satellites, 
e.g. comets and asteroids. If the population of these objects 
is similar to that around the Sun, there must be a consider- 
able number of soli d objects tha t mov e around the GC. In 
our previous paper (jCadez et al.l (|2008l ). hereafter CCK08), 
it has been shown that tidal evolution of orbits can bring 
their periastron close to the black hole. If an object finds 
itself below its Roche radius, tidal forces dominate and dis- 
rupt it. In this paper, we study the possibility that tidal 
disruption events are connected with the observed flares 
from the GC. The timescales of flares strongly suggest mo- 
tion only a few Schwarzschild radii away from the black hole 
and it has been speculated that the energy released during 
these flar es corresponds to t he source mass of the order of 
~ 10 20 g (|Genzel et al.ll2003t) . Asteroid-like objects explain 
both observational results: they are of the correct mass and 
are tidally disrupted closer to the black hole than gaseous 
blobs. In Active Galactic Nuclei, these events would be un- 
detectable, but in the inactive centre of our Galaxy, it is 
possible to observe them. 
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In Sect. [21 we review some ideas regarding the popula- 
tion of these small bodies in the immediate neighbourhood 
of the Galactic centre black hole and discuss the difference 
in the Roche radius for solid and gravity-dominated objects. 
In Sect.[31 we study tidal disruption of a small object nu- 
merically, using ray-tracing techniques in a Schwarzschild 
background. We present calculations of the temporal evo- 
lution of the light curve of infalling objects as a function of 
various parameters. In Sect. [4] the possible consequences of 
tidal squeezing on the object's internal magnetic field are 
discussed, and calculated light curves are compared with 
the observed light curves of flares in the Galactic centre. 

2. Small bodies at the Galactic center 

It appears reasonable to assume that stars at the GC have 
planets and other small orbiting bodies, such as asteroids 
and comets. When the parent star approaches the central 
black hole, tidal interaction may either strip these bodies 
off their parent star or cause them to become more bound. 
In this way, the star may lose its satellites before it is swal- 
lowed or disrupted by the black hole. Little is known about 
the dynamics that determine the fate of these satellites. 

A related problem was studied by iBerukoff fc Hansen] 
(20061), who investigated cluster core dynamics in the GC, 
by considering intermediate massive black holes that trans- 
port multiple tightly-bound stars inwards into the potential 
well to a tidal-stripping radius. They found that this mech- 
anism is likely to produce a population of stars close to the 
GC on high-eccentricity, high- inclination orbits. An analo- 
gous mechanism could be applicable to a planetary systems 
of stars. 

The possibility of detecting gravitational waves with 
LISA stimulated an extensive investigation of mechanisms 
to c ause stars to move inwards into a massive black 
hole (|Hopmanll2006l and references therein). These studies 
demonstrated that the central density cusp of the stellar 
distribution flattens to n(r) oc r -1 / 2 at a critical radius of 
oqw = few x 0.01 pc for the distribution of 1 M Q objects. 
Here oqw was defined to be the radius where gravitational 
radiation damping prevails over resonant relaxation in the 
central star cluster. Since the power of gravitational wave 
emission is proportional to to 2 , acw is negligible for small 
bodies. Consequently, their central density cusp may ex- 
tend much further to radii at which energy and/or angular 
momentum can be extracted efficiently from the orbit by 
tidal interaction, as discussed by CCK08. 

A beautiful presentation of the Galactic centre envi- 
ronment based on hig h resolution observations is given by 
iPaumard et a l. (2006). They confirm a marginally steeper 
stellar density cusp n(r) oc r -2 , and its flattening, in fact a 
sharp drop, at ks 0.04 pc. Yet the central cusp is not located 
in a spherically symmetric configuration but in two young 
stellar disks within 0.5 pc of Sgr A* containing a mass of 
- 1.5 x 10 4 M Q . 

Further clues about the causes and effects of tidal evo- 
lution in a many-body system may come from detailed dy- 
namical studies of the evolution of the Solar system, which 
have revealed a number of unexpected phenomena that 
perturb orbits and tidally melt smal ler bodies when they 
are caught in stochastic resonances (iDermott et aTJll988t 
iFaber et al.ll2005t iMudrvk fc Wull2006fh 

As for the poss ible number of t hese satellites, we note 
that according to ICochran et al.1 (|1995l ) the Edgeworth- 



Kuiper belt of our Solar system may still (4.5 billion years 
after its creation) contain as many as 2 x 10 8 objects of 
radii < 10 km. 

On the basis of this discussion, one might expect to find 
stellar-system satellites all the way down to the black hole, a 
fair proportion of them on low-periastron, highly-eccentric 
orbits. 

Close to low periastron, tidal forces apply significant 
work to objects. This work partially dissipates (in terms 
of heat, or accelerated electrons), lowering the orbital en- 
ergy and initiating significant tidal evolution of the orbit, 
CCK08. Tidal interaction is strongest in resonance, i.e. 
when the orbital and fundamental quadrupolc frequencies 
are identical. In this context, one must consider two classes 
of these satellites: those that are gravity-dominated and 
those that are solid-state-forces-dominated. For gravity- 
dominated objects, the fundamental frequency is u g w 
2y/ Gp/3n, while for solid-state-dominated objects, it is 
v s = jc s /R (p is the density of the body, R its radius 
and c s the speed of sound). If we assume typical values 
of c s ~ 5 km/s and p w 5 g cm -3 , we ascertain the ra- 
dius dividing the two classes to be about 1000 km, i.e. the 
radius of the asteroid Ceres. Thus, all gravity-dominated 
satellites will have about the same fundamental quadrupolc 
frequency, corresponding to the period of about 54 min- 
utes. All smaller satellites will have shorter fundamen- 
tal periods. Gravity-dominated satellites begin their rapid 
tidal evolution when their periastron reaches r p w 10 r g 
(r g = GMbh/c 2 ), while solid-state-dominated bodies may 
begin significant tidal evolution even closer to the black 
hole. 

To illustrate the difference between the two classes 
of objects, we introduce the Roche radius and the ef- 
fective Roche radius. The effective Roche radius is de- 
fined as the radius at which uj q T{ = 1, where w q is 
the fundamental resonant frequency of quadrupole modes 
and Tf = ^r^/GM^) 1 / 2 is the periastron crossing time 

(|Gomboc fc Cadezl 120051 ) Fl FigureQ] shows the Roche ra- 
dius and the effective Roche radius as a function of object 
size for different masses of gravity-dominated and solid- 
state-dominated objects. The positions of different celestial 
objects are also marked in this diagram. We note that the 
effective Roche radius for solid-state-dominated asteroids 
is below 2 r g , while the Roche radius for these objects is 
about 10 r g . 

3. Tidal disruption of a small body and resulting 
light curves 

We are interested in small objects orbiting the black hole 
on eccentric, low-periastron orbits, which are presumably 
brought there by processes described in Sect. [2] An object 
orbiting a black hole on such an orbit experiences high 
tides and is therefore heated at each periastron passage. 
As we showed in CCK08, for a range of object masses, 
solid-state forces, gravity, and heat conduction interact in 
an interesting way. Namely, the varying stress induced by 
the varying tidal field exerts work. This work is transferred 
into heat that escapes on a thermal diffusion timescale 
r solid = R 2 tt~ 2 D~ 1 , of the order of 3 x 10 5 years for 10 km 



The effective Roche radius becomes the Roche radius for 
gravity-dominated objects. 
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Fig. 1. Roche radius and effective Roche radius as a func- 
tion of mass and radius of the object. Blue lines indicate 
the Roche radius for gravity-dominated objects and red 
lines the effective Roche radius for solid-state-dominated 
objects (Mce is the mass of the asteroid Ceres). Lines of 
constant sound velocity c s are plotted in light grey and 
refer to c s = 1.5, 3.5, 5 km s _1 (ice, rock, iron). Some 
celestial objects are marked by blue dots: Earth, Jupiter, 
Sun, Sirius B, Aldebaran, Betelgeuse and Rigel. The green 
line represents zero-age main-sequence stars and the purple 
line represents white dwarfs. 



objects consisting of normal, solid rock (thermal diffusion 
D ~ 10 _6 m 2 s~ 1 ), which is long with respect to the orbital 
period. Thus, a stress that is well below breaking stress 
can do enough accumulated work during the thermal diffu- 
sion time to heat the body to melting temperatures. When 
melting begins, the solid-state structure weakens, allowing 
tides to do more work per passage, which leads to complete 
melting in a short time afterwards. In CCK08, we assumed 
that an object is on a sufficiently eccentric orbit (e > 0.1) 
and that its fundamental resonance has the mechanical Q 
value of 500. We found that objects with R > 20 km melt 



if r p < 5 r e 



and those with R > 50 km melt if r p < 10 r B 



However, heat produced by tides is insufficient to evaporate 
small bodies, because the weakening of solid-state structure 
allows convection to set in so that heat is efficiently brought 
to the surface, where it is radiated away at sub-evaporating 
temperatures. In this way, small solid-state-dominated ob- 
jects that were orbiting the black hole above their effective 
Roche radius, find themselves below their effective Roche 
radius, simply because melting has moved their effective 
Roche radius closer to the Roche radius. This transition 
occurs without appreciable loss of orbital energy. It will 
be shown below that evaporation of the body occurs only 
during the last runs when the orbit develops into into an 
unstable, critical orbit. 

In previous simulations, hot spots orbiting a 
black hole were treated either as point sources 
(iKaras fc Bad Il992t Hao et all Il994t ISchnittmanl 12004 



Pechacek et al. 20061) or small, fixed- size blobs (|Karasl 



1996s iSchnitt man fc Bertschingerl 12004}) . However, as we 



argued above, even solid objects larger than 10 or a few 
10 km melt when their periastron is below 10 r g , and 
therefore rapidly deform afterwards. 



Here, we study the final stages of the tidal capture of 
such a small, melted object orbiting a Schwarzschild black 
hole. 

To ensure that the calculations are manageable, we 
study the following simplified problem. We consider a solid, 
spherical object moving around a black hole on an elon- 
gated, bound orbit. After tidal forces do sufficient work, 
the object momentarily melts. The distance from the black 
hole at which this happens is denoted by r mo i t . After this 
moment, pressure forces become negligible with respect to 
inertial forces, and the object begins to behave as a collec- 
tion of freely noninteracting particles. We use ray-tracing 
from points on the object's surface to the distant observer 
to determine its appearan ce and light curve. Ray -tracing 
procedures are described in lCadez fe Kosticl (|2005l) and are 
freely available from the authors. 



3.1. Numerical code 

The surface of the object is triangulated with respect to 
the mesh of spherical coordinates on the original sphere. 
Before melting, the sphere is considered to be rigid and 
non-rotating (zero-spin angular-momentum). Its centre of 
mass follows a time-like geodesic with a chosen specific en- 
ergy E/m and specific angular momentum l/m. The coor- 
dinates of mesh points are expressed in the Fermi- Walker 
transported frame local to the centre of mass and are ini- 
tially constant with respect to proper time. The rigidity 
constraint prescribes the same velocity for each point on 
the sphere and thus different specific energy and angular 
momentum. When the sphere is melted, it finds itself be- 
low the Roche radius. From this moment on, the constraint 
no longer applies and surface mesh points begin to follow 
time-like geodesies specified with their inherited specific en- 
ergy and angular momentum. 

The object is considered as a monochromatic radiator 
(frequency constant rest-frame emissivity Iq) emitting 
according to Lambert's law in its own rest frame. The image 
of the object is constructed by mapping triangles on the sur- 
face of the object onto the image plane of the observer and 
assigning the appropriate intensity and frequency to their 
image. As long as triangles on the object are sufficiently 
small and gravitational lensing is moderate, their images 
are also triangles that are defined by the positions of their 
vertices (see Fig. [5]). We note, however, that as time goes on, 
the stretching of the object increases exponentially, so that 
sooner or later every triangulation becomes too scarce. In 
such cases, we halt the simulation. Through gravitational 
lensing, each vertex on the object produces multiple images. 
We note that the size of higher order images exponentially 
decreases with order, so it is sufficient t o take into account 
only the primary and secondary image (jCadez fc Gombocl 
fl99l . 

The image is observed in the image plane at constant 
time intervals with respect to the observer. Since we are 
dealing with a highly dynamical problem, we must take 
into account that light travel times from the object to the 
observer are different for rays from different vertices. To 
calculate the position of the image of the vertex say Vi at 
observer time i bs, we first calculate the trail of the image of 
this vertex. We define Ci = {r, rj(r), 9%(r), 0i(t)} to be the 
orbit of the point Vi, parametrized with Schwarzschild co- 
ordinate time r. We send rays from Vi toward the observer 
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Fig. 2. Mapping the triangulated sphere onto its image. 




Fig. 3. Determining the trails and c\ 2 ^ of the image of 
a vertex Vi and its position at t Q h s and i Q bs + At. 



and calculate the trails C- = {t[^ (r), (r), y t - (r)} and 
Cj- 2 ^ = {t- 2 ^(r),a;- 2 ^(T),?/- 2 ^(T)} of the primary and sec- 
ondary image. Here t\ a \r) = t + At\ a ' and At^ is the 
time of flight from Vi to the observer along the a-th ray 
(a = 1,2). The orbit is calculated at discrete time inter- 
vals T k and the coordinates in the image {xi,yi} at time 
fobs are calculated by linear interpolation as indicated in 

Fig. [3] The intensities /■ , l[ 2 ^ and frequencies v % , 
are interpolated in a similar way. In the image, all vis- 
ible triangle areas Si are calculated and they contribute 
AL;(f bs) = h (fobs) Si (fobs) to the total observed luminos- 
ity. 

3.2. Results 

We simulated objects moving along different orbits, charac- 
terized by different values of E/m and l/m. We found that 
relativistic tidal effects depend predominantly on the ratio 

C = (E/mc 2 - V min )/(V max - K min ) , (1) 

Where V max and V m i n are the maximum and min- 
imum of the relativistic effective potential V = 

((l-2r g /r)(l + (l/mcr) 2 )) 1/2 dMisner et alJ 11973ft . For 
stable circular orbits, £ = 0. Tidal effects on these orbits 
are small and evolve as in the Keplerian case, where tidal 
deformations (below Roche limit) increase linearly. On the 
other hand, £ ~ 1 corresponds to semi-stable orbital that 

2 Such orbits partly wind around the black hole at radii in the 
interval 3r g < r < 6r g , make short excursions away from the 
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Fig. 4. Area S of the image of a tidally-evolving object seen 
above the orbital plane as a function of time for different or- 
bital parameters: a) elliptic orbit close to innermost stable 
circular orbit (ISCO) with r A = 6.2122 r g , r P = 6.2097 r g 
and ( = 2.44 x 10 -5 ; b) plunging parabolic orbit with 
E/mc 2 = 1 and l/mr g = 3.999998 and ( = 1 + 6.6 x 1CT 6 ; 
c) plunging hyperbolic orbit with r(V max ) — 3.6 r g and £ = 
1 + 1(T 6 ; d) elliptic orbit close to ISCO with r A = 6.1 r g , 
rp = 5.95 r g and ( = 1 — 7 x 10~ 6 ; e) elliptic orbit with 
r A = 20 r g , r P = 4.45 r g and C = 1 - 7 x 10~ 6 ; f) scat- 
tering parabolic orbit with E/mc 2 = 1, l/mr g — 4.000002 
and C — 1 — 7 x 10~ 6 ; g) plunging hyperbolic orbit with 
E/mc 2 — 3.0, l/mr g = 6, £ = 6.35; h) linearly growing 
deformation. Sq is the visible area of originally rigid sphere 
as seen by the observer. In cases b,c,e and f, the simulation 
was halted before the object deformed beyond a limit of 
numerical accuracy. 

extend to Vmax- Simulations exhibit exponentially growing 
tidal deformations (see Figs. 2] and [5]), so that even small 
bodies, which generally suffer less, are eventually exponen- 
tially smeared into long, thin threads. The extent of tidal 
deformation is estimated from the length of the object's 
image as seen by the observer above the orbital plane. In 
our simulations, all the exponential timescales become sim- 
ilar, r c « 11 r g /c. To assess the dependence of the speed 
of tidal deformations on the initial melting radius r me it, 
we simulated motions along the ( = (1 — 7 x 10~ 6 ) orbit. 
As shown in Fig. [51 tidal effects do not depend strongly on 
fmcit- We note that the Keplerian case is quite different: if 
r mc it is close to periastron (below the Roche radius), tidal 
deformations increase linearly with time, while they do not 
increase at all if r me it is close to apastron (above the Roche 
radius). 

Here, we prese nt light curves for two limiting cases: £ ~ 
1 and £ -> (see iKostid (|2008ft for a detailed discussion . 
Time initiates from the moment of first detection by the 
observer. 

The £ ~ 1 case is represented by a parabolic orbilQ 
with E/mc 2 = 1 and l/mr g c = 3.999998 (case 'b' in Fig. 
[4j. The starting point is at ro = 50 r g and r me it = r o- 

black hole, and repeat the winding cycle or end up behind the 
horizon. They do not exist in the Keplerian case. 

3 A parabolic orbit in this context is meant as a synonym for 
highly eccentric elliptical orbit, while motion along hyperbolic 
orbits is considered as a test of proper behaviour of numerical 
models. 
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Fig. 5. The length / of the image of a tidally-evolving ob- 
ject seen above the orbital plane as a function of time for 
different orbital parameters: a) elliptic orbit close to in- 
nermost stable circular orbit (ISCO) with ta = 6.2122 r g , 
r P = 6.2097 r g and ( = 2.44 x 10~ 5 ; b) plunging parabolic 
orbit with E/mc 2 = 1, l/mr g = 3.999998 and ( = 1 + 1CT 6 ; 
c) plunging hyperbolic orbit with r(V max ) = 3.6 r g and ( — 
1 + 1CT 6 ; d) elliptic orbit close to ISCO with r A = 6.1 r g , 
rp = 5.95 r g and ( = 1 - 3 x 10 -6 ; e) elliptic orbit with 
r A = 20 r g , r P = 4.45 r g and ( = 1 - 3 x 10" 6 ; f) scat- 
tering parabolic orbit with E/mc 2 = 1, l/mr g — 4.000002 
and C = 1 — 3 x 10 -6 ; g) plunging hyperbolic orbit with 
E/mc 2 — 3.0 and l/mr g — 6 and ( = 6.35; h) linearly 
growing deformation. In cases b,c,e and f, the simulation 
was halted before the object evolved beyond a limit of nu- 
merical accuracy. 
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Fig. 7. Top: Orbit of the object with E/mc 2 — 1 and 
l/mr g c = 3.999998. Bottom: The light curve of the ob- 
ject of size R = 1Q~ 5 r g . The "melting point" is at 
r mcit = r = 50 r g . The observer is 1° above the orbital 
plane. Lq is the luminosity of originally rigid sphere. 
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Fig. 6. Area S of the image of a tidally evolving object seen 
above the orbital plane as a function of time for £ = 1 — 7 x 
10~ 6 parabolic orbit (case 'f from Fig. [4]) with r$ — 100 r g 
for different r mc i t . So is the visible area of originally rigid 
sphere as seen by the observer. Note that the dependence 
on r me it/V g is weak 



Two examples of light curves are shown in Figs. [7] and [5] 
with corresponding images of the object in Figs. l8l and [T0l 
By comparing the top right images of Figs. 151 and [T0l one 
realizes that in the case of Fig. [51 where the observer is 
only 1° above the orbital plane, the small object has been 
gravitationally lensed into primary and secondary Einstein 
arcs. This magnification has a duration of about 1 r g /c. In 



the last two images of Fig. the object is already stretched 
into a thin thread as confirmed in Fig. 1 101 

The sharp and narrow peaks in the light curves are due 
to gravitational lensing and always appear in pairs: the 
lower one is the contribution of the secondary image and 
precedes the one from the primary image. The peaks are 
higher if the observer is closer to the orbital plane (compare 
with Fig. [5]). The wider bumps following sharp peaks are the 
result of Doppler boosting and the aberration beaming of 
light. The overall increase in luminosity is a consequence of 
the exponentially increasing surface of the object. 

To compare results with orbiting blobs of fixed size, 
we also show light curves of a freely falling rigi d object in 
Figs.F Fland p The s e resu lts are in agreement with lBao et al.l 
(|1994[ ) and lKarasl (|1996| ). Both light curves include the sec- 
ondary image contribution. 

The case C ~ * is illustrated by tidal evolution along 
a quasi-circular orbit with E/mc 2 = 0.94294053 and 
l/mr g c — 3.46610164, so that the apastron lies at r A — 
6.216 r g and periastron at rp ~ 6.206 r g (case 'a' in Fig. 
U). The starting point lies at = 6.21 r g and r me i t = r A . 
The corresponding light curves for different observer's ori- 
entations are shown in Fig.[TT] In Fig.rTJl we also show 
the light curve of the same object after a large number 
of periods. Since the light curve is not smeared out and 
extremely sharp spikes are still present at the end, we con- 
clude that the object is not completely tidally disrupted; its 
size along the orbit increases (linearly) from R = 10~ 5 r g 
to e.g. R ~ 0.1 r g . 
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TIME=376.000 r g l c |TIME=402.000 r g l c 



Fig. 8. The image of a freely-falling surface that was ini- 
tially spherical. It is seen by an observer 1° above the orbital 
plane. The colours on this tidally-distorted, gravitationally- 
lensed image correspond to redshift. The green circle rep- 
resents the Schwarzschild radius. 



Fig. 10. Freely-falling spherical surface, tidally-distorted 
seen through gravitational lensing, as seen by the observer 
70° above the orbital plane. An enlargement of the pri- 
mary image is included in the top left image. The colours 
correspond to the redshift. The green circle represents the 
Schwarzschild radius. 
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Fig. 9. Top: Orbit of the object with E/mc 2 = 1 and 
l/mr g c = 3.999998. Bottom: The light curve of the ob- 
ject of size R — 1CP 5 r g . The "melting point" is at 
'"melt = r a = 50 r g . The observer is 70° above the orbital 
plane. Lq is the luminosity of originally rigid sphere. 
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Fig. 11. Top: Orbit of the object with E/mc 2 = 0.94294053 
and l/mr g c — 3.46610164. Bottom: Light curves of the ob- 
ject of size R ~ 10~ 5 r g for different angles 9 of the observer. 
Lq is the luminosity of originally rigid sphere. 
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Fig. 12. Same as Fig.Ql] for 9 = 89° of the observer but 
only at a much later time. 

4. Discussion 

4.1. The effects of magnetic field 

Most authors now agree that radio and NIR flaring emission 
from Sgr A* is produced by synchrotron radiation of elec- 
trons with Lorentz fact ors 7 e ~ 10 2 — 10 3 moving in a mag- 
netic f ield (10-100 G) (fMarkoff et al] l200HlBaganoff et all 
l200lt I Yuan et al.l |2004 ). while X-rav flares are assumed 
to be produced by inverse Compton sc attering of low en- 
ergy p hotons from relativistic electrons (Eck art et al.l [2004. 
2006a). As to the observed magnetic field, there is still 
no general agreement on its strength at the center of the 
Galax y. However, it has been pointed out bv lLaRosa et alJ 
(2006J) that the dispute is whether the field is strong (of 



order 1 mG) and globally organized (jMorris &; Serabvnl 
Il996t iMorris! l200ll . or globally wea k (of order 10 uG) 
with regions in which it is st ronger (jLaRosa et alJ l2005t 



iBoldvrev fc Yusef-Zadehl ^006). In any case the observed 
strength is orders of magnitude lower than requested by 
synchrotron models. Assuming that the object is perme- 
ated with such a low magnetic field, our simulations sug- 
gest a simple mechanism for increasing its strength as well 
as providing electrons at high j e . 

Our simulations, discussed in the previous section, 
demonstrate that objects that find themselves on ( ~ 1 
orbits become exponentially squeezed and elongated. Since 
the inward squeezing velocity is low (v/c ~ 10 -7 ), it is 
possible to use the MHD approximation for describing the 
evolution of the magnetic field 



-Vx(VxB) = -f 



V x (v x B) , (2) 



where a is the electrical conductivity and v is the plasma 
velocity. With typ ical values a = 4.14 x 10 6 Q _1 m _1 for 
molten aluminum (jRhim fc Ishikaw ail998l ) and L — 100 km 
for the size of the object, the magnetic Reynold's number 
R m becomes large 

Rm = \Iq<JvL Ri 10 . 

Consequently, the diffusion term of Eq.[5] becomes negligi- 
ble, and we can write 



which, according to Alfven's frozen-flux theorem, leads to 
magnetic-flux conservation. As the object with initial mag- 
netic field Bq and cross section Sq is exponentially squeezed, 
the magnetic field increases likewise as B = B So/S. 
Taking i nto account that t he tidal deformation tensor is 
traceless dMisner et al.ll973l) , the average cross-section area 
of a squeezed and elongated object is inversely proportional 
to its length, so that an approximate magnetic field at a 
later time can be estimated to be 



B(t) » B j^-Av(t) 
4 Ho 



(4) 



where r is the distance from the black hole, Ro the initial 
size of the object, and Aip the arc length of the object. The 
exponentially increasing value of Atp(t) can be either calcu- 
lated from tidal acceleration (Sect. I3T2")) or inferred directly 
from pictures, after taking into account relativistic length 
corrections. 

At the beginning, we assume that the object is perme- 
ated with a (galactic) magnetic field Bq = 1 mG, has a 
size of Rq — 50 km w 10~ 5 r g , later orbits around the 
black hole at r = 4 r g and is stretched into an arc with 
Aip — 10°. After inserting these numbers into Eq.QJ the 
estimated magnetic field can become as high as B w 50 G. 
We also note that the initial magnetic field of the object 
may well be far higher than the assumed value of 1 mG, 
since magnetic fields as high as Earth's magnetic field have 
been measured in asteroids ([Kivelson et al.l H993). 

We also highlight that the exponentially-increasing 
magnetic field satisfies the conditions within a betatron. 
Therefore, by taking into account the initial magnetic field 
and the magnetic field rise-time, it is not difficult to obtain 
electrons of extremely high 7. 



4.2. Flares 

We consider whether the flare activity observed at the GC 
is due to the capture of a small body, of comparable mass to 
that of a comet or an asteroid. The only important physical 
constraints on the source are the requirement that: 1) the 
source can provide sufficient energy, and 2) it is sufficiently 
small to produce observable modulation of the light curve. 

The timescales in our simulations agree with the 
timescales of flares, if we assume that the source of the flares 
moves on an unstable, circular orbit at r < 4 r g around a 
Schwarzschilcfl black hole of mass 3.6 x 10 6 Af©. Since this 
is a £ ~ 1 orbit, the tidal work increases exponentially 
and the object is deformed exponentially on a timescale of 
11.3 r g /c. In this case the timescale for magnetic diffusion 
is much larger than the dynamical timescale, therefore the 
magnetic flux is conserved. Consequently, as the object is 
squeezed, the magnetic field (exponentially) increases and 
the electrons with high 7 begin to emit synchrotron radia- 
tion. 

We assume that the energy released during this process 
originates in the gravitational potential energy of the object 
and is thus proportional to its mass. Since potential energy 
differences in an orbit at r < 4 r g are of the order of a 
few percent of mc 2 , it follows that objects producing them 



dB 

~dt 



V x (v x B) 



(3) 



4 All analyses has been completed in Schwarzshild space-time, 
but in order to fit light-curves it is sometimes preferable to con- 
sider Kerr space-time with a moderate value of the angular mo- 
mentum parameter a. 
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proba bly have masses of the order of 10 20 g (|Genzel et al.l 
2003). If the sources of flares were gaseous blobs of such 
a large mass, they would find themselves below the Roche 
radius far away from the black hole, and would therefore 
be completely disrupted before producing any modulation 
of the light curve. We conclude that the source of the flares 
is a small, solid object that orbits the black hole above the 
effective Roche radius. 

After many orbits, this object melts due to work exerted 
by tidal forces, thus its effective Roche radius increases 
quickly towards the Roche radius. As the melted object 
approaches the black hole on a £ ~ 1 orbit, it is deformed 
into a long, exponentially increasing, tidal tail (similar to 
tails in Figs.l8l and [TOf as it falls within the Roche radius. 
Since the size of the object is initially small, it can still pro- 
duce the modulation of the signal for a number of orbits, 
until it is completely smeared along the orbit. 

Using ray-tracing techniques, the appearance and the 
luminosity of such an event as seen by observers was cal- 
culated at different inclinations and position angles with 
respect to the orbit of the infalling asteroid. In this sim- 
ple toy-model, we introduced five parameters: mass of the 
black hole Mbh, tidal heating timescale Th, length of the 
tidal tail /tail, inclination of the orbit i, and the longitude 
of the line of nodes f2. A robust description of the June 
16th 2003 IR flare was found for the following parameter 
values: M bh = 4 x 1O 6 M0 , r h = 2300 s, l taa = 10 8 km, 
i « 90° (jCadez et al.l l2006h . Despite the model simplicity, 
the value of Mbh that we find is in good agreement with 
the measured value. 

We were also able to describe the XMM-Newton /EPIC 
light curve of the flare of April 4, 2007 (jPorquet et al.ll2008l ) 
with the following scenario: after the tidal deformation de- 
velops, the deformed object can be imagined as a long string 
of freely falling beads threading the same orbit. The light 
curve S(t) of such a string can thus be described by a con- 
volution 



S(t) = / G(t,t')K(t')dt' 



(5) 



where G(t, t') is the Green's function belonging to the orbit, 
i.e. it is the signal, as a function of time t, of a single bead 
crossing a fixed reference point at time t', and K (t') is the 
distribution of luminosity among the beads. The Green's 
function depends on the orbital angular momentum I and 
energy E of the orbit. We are interested in plunging orbits 
with C ~ 1) since tidal evolution produces these orbits, and 
only in such orbits can the length of the tidal tail grow 
exponentially and produce a flare. For such an orbit, the 
Green's function is dominated by a series of sharp spikes 
produced during the last few turns before plunging behind 
the horizon (see Fig. [7]). The quasi-period of the spikes de- 
pends on the position of V max . We use the observed signal 
S(t) and solve Eq.[5]for the luminosity distribution K{t'). 
Depending on the Green's function used, different solutions 
are obtained, most of which are physically unrealistic highly 
variable functions. Only for l/mr g c = 4.3 and E — 1.04 mc 2 
does the distribution become a simple positive function 
K(t'), which can be well approximated with a steep, expo- 
nential rise, followed by a slower, exponential decayO Both 
K(t') and its smooth approximation K{t') are shown in the 
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5 The energy and angular momentum were calculated by as- 
suming that Mbh = 3.6 x 1O 6 M . 



Time [h] 

Fig. 13. The X MM-Newton/EPIC li ght curve of the flare 
of April 4, 2007 (jPorquet et al.ll2008D . fitted with the sim- 
ple intensity distribution shown in the inset at top right. 
Residuals are shown as crosses at the bottom. 



inset of Fig.[T5J The function K{t') and the Green's func- 
tion described above thus represent a viable description of 
a flare produced by an originally small object that is tidally 
distorted into a long thread during the last exponential part 
of its evolution in a ( » 1 plunging orbit. The exponential 
rise time of the flare is 4.8 min and the exponential decay 
time is 12.5 min, which translates into a tidal-tail length of 
^taii = 1.1 x 10 8 km. The heating time that was estimated 
in the case of the IR flare cannot be estimated in this case, 
because we only have two peaks in the Green's function. 
In the case of the IR flare, the exponential rise time was 
considered as being rapid, and, considering the noise in the 
data, to be consistent with the short rise time preferred by 
the X-ray fit. The orbital energy obtained is probably un- 
reasonably high. We note, however, that our analysis was 
performed in the Schwarzschild space-time. If this analy- 
sis were completed instead in Kerr space-time, the Green's 
function quasi-period would decrease for co-rotating orbits. 
For example, if the value of the Kerr parameter was 0.4, the 
orbital energy would be only w 0.945 mc 2 . 

On the other hand, if the object orbits the black hole on 
a ( ~ orbit (i.e. ISCO), tidal work and deformation in- 
crease only linearly. In this case, the timescale for magnetic 
diffusion is much shorter than the dynamical timescale, the 
magnetic flux is therefore not conserved and the magnetic 
field decays before reaching sufficiently high values to pro- 
duce synchrotron radiation. 

5. Conclusions 

A puzzling characteristic of the Galactic centre that har- 
bours a massive black hole are Galactic flares detected 
at various wavelengths. Their main characteristics are: a) 
short duration, b) quasiperiodic oscillations of period 17-22 
minutes, c) total energy release of the order 10 39 - 5 erg, and 
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Fig. 14. What the XMM-Newton/EPIC flare of April 4, 
2007 might look like, if it could be seen from above the 
orbital plane. 



d) strong indications of the presence of magnetic fields of 
up to 100 gauss. 

The exact cause of the flares is still unclear, and several 
models have been proposed to explain their origin (see ref- 
erences in introduction). In this paper, we have proposed 
and explored the idea that flares are produced by the final 
accretion of dense objects with mass of the order of 10 20 g 
onto the massive black hole in SgrA*. We have studied the 
effects of black- hole tides on small, solid objects in the vicin- 
ity of a massive black hole. Such objects are expected to 
populate the GC as the result of their being stripped from 
their parent stars. In CCK08, we showed that solid objects 
in the mass range ~ 10 19 — 10 21 g melt in the near vicinity of 
the black hole. Afterwards their orbital evolution naturally 
leads to capture orbits, most likely by means of unstable, 
circular orbits. In an exhaustive numerical simulation, we 
studied the evolution of shapes of objects and their light 
curves as the objects move toward the event horizon of the 
black hole. We found that tides generally stretch the ob- 
jects into long, thin threads that extend exponentially on 
a timescale r c ~ 11.3 r g /c, if the orbital criticality param- 
eter £ is close to 1. Furthermore, if the object is larger 
than > 0.3 r g , it is completely disrupted before producing 
any significant modulation of the light curve. Therefore, the 
object should be relatively small to produce the observed 
modulation. 

Tidal evolution of the shapes and the images of these 
objects was then studied numerically in the space-time of 
a Schwarzschild black hole. Because of the complexity of 
the problem, we did not consider the more general Kerr 
space-time and instead limited the analysis to pressureless 
gas. We found in general, that these objects are eventually 
elongated into thin threads that in £ ~ 1 orbits extend 
exponentially. 

During this exponential elongation, we highlighted that 
conditions for magnetic-flux conservation were obeyed, so 
that relatively high magnetic-field strengths could be gener- 
ated. The exponential growth of the magnetic-field density 
also allows the betatron mechanism to operate, and possi- 
bly generate highly relativistic electrons. We propose this 
process as a possible source of Galactic flares. We have not 
yet developed the precise mechanism for extracting radia- 
tion from gravitational potential energy, but we are able to 



model light curves of Galactic flares with simple assump- 
tions. 

Finally, we emphasize that our scenario differs from 
other models in the following respects: 1) Flares origi- 
nate in marginally-bound and not marginally-stable or- 
bits. This allows them to travel closer to the black hole 
and provides a mechanism for the rapid extraction of en- 
ergy. It also agrees with measurements of the size of the 
flaring region in the Gala ctic centre, which imply that 
it is not larger than 6 r g (|Reidl [2008). and relaxes con- 
straints on the angular momentum of the black hole. 2) 
Our model does n ot require a fa i rly ra pi dly rotating black 
hole a s nee ded by I Genzel et all (|2003f ). lAschenbach et all 
(|2004l ) and iTrippe et all (j2007f ). 3) An accretion disk at 
the G alactic centre is not required. lYusef-Zadeh et all 
(|2008f ) highlited that "there is no definitive evidence that 
Sgr A* has a disk" . 4) Elongated structures inv o ked i n 
other models to fit flare data (e.g. iFalanga et ail (120071) ; 
IZamaninasab et all (|2008l ) and iHamaus et all (|2008l ))"~de- 



velop naturally in the proposed model as a result of the 
tidal evolution of a melted body (see Fig fLl")) . 5) The evolv- 
ing elongated structures also provide a natural explanation 
of rapid changes in magnetic-field density inferred from ob- 
servations, increasing first as the body is squeezed by tides 
and finally decreasing as the body crosses the horizon of 
the black hole. 
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